Abstract. We have analyzed the ultraviolet and optical 
emission line ratios of a large sample of extragalactic 
radio sources (QSOs and radio galaxies), with the help 
of models combining AGN photoionization and shocks. 
The results strongly suggest that the two ionizing 
mechanisms frequently coexist. The model sequences 
obtained by varying the balance between shocks and 
AGN photoionization account for most emission line 
data in the 12 line ratio diagrams we have considered. In 
the frequently used diagrams involving [OIII]A5007/H/3, 

^ [OI]A6300/Ha, [NII]A6584/Ha and [SII]AA6716,6731 
■ (Veilleux & Osterbrock 1987), the effect of varying the 
shock-photoionization balance mimics a variation of the 

(N ionization parameter (U) in traditional photoionization 



> 



sequences. In most of the remaining diagrams, such as 
[OI]A6300/[OIII]A5007 vs. [OIII]A4363/[OIII]A5007, 
[OIII]A5007/H/3 vs. [OIII]A4363/[OIII]A5007 and 
CIII]A1909/[CII]A2326 vs. CIVA1549/[CII]A2326, the 
data can only be accounted for if both photoionization 
and shocks contribute to the line fluxes. The coexis- 
tence of shocks and AGN photoionization also provides 
an explanation for the most extreme objects in the 
NVA1240/HeIIA1640 vs. NVA1240/CIVA1549 diagram 
. without requiring largely super-solar metallicities. In 
j addition, we show that there is a relationship between the 
O ■ [OII]A3727/[OIII]A5007 ratio (i.e., the ionization level of 
' the gas) and the radio size in radio galaxies. This strongly 
O . supports the hypothesis that the most compact (< 2kpc) 
I ' and the largest (> 150 kpc) sources are dominated by 
^~ . photoionization, while intermediate-sized radio galaxies 
' are dominated by shocks. We briefly discuss the possible 
C$ • origin of the relation between the shock- ionization balance 
^ ' and the radio size. 

•i-H . 
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1. Introduction 

One of the most debated questions about the interaction 
between active galactic nuclei (AGN) and their environ- 
ment is the nature of the ionizing source of the extended 
structures of ionized gas detected in powerful radio galax- 
ies (Spinrad & Djorgovski 1984; Tadhunter et al. 1988; 
Baum et al. 1988) and QSOs (Hes et al. 1996). The emis- 
sion line spectra of these regions impose the strongest con- 
straints on this issue. Initial studies of extended line emis- 
sion in radio sources showed that photoionization models 
provide an acceptable fit of emission line ratios in the op- 
tical (Robinson et al. 1987; Saunders et al. 1989; Prieto 
et al. 1993) and the UV (Villar-Martin et al. 1997), in 
good agreement with the "unified scheme" (e.g. Antonucci 
1993). 

However, some properties of the extended emission 
cannot be accounted for by the photoionization hypoth- 
esis. Some AGN have a spectrum dominated by low ion- 
ization lines (e.g. Koski & Osterbrock 1976; Durret 1990), 
or show [OIII]A4363/[OIII]A5007 ratios indicative of high 
electron temperatures (T c ~ 20000 K; Ferland & Oster- 
brock 1987). These two properties are more reminiscent 
of ionizing shocks. Shocks can be produced, for example, 
by the interactions of the radio ejecta with the interstel- 
lar medium (ISM), as suggested by the spatial coincidence 
between radio and line emission in radio sources (Baum & 
Heckman 1989; McCarthy et al. 1995; Ridgway & Stock- 
ton 1997; Lchnert et al. 1999). The disturbed kinemat- 
ics of radio galaxies are also indicative of the influence of 
shocks in the ionizing process (Villar-Martin et al. 1999a). 
With the help of theoretical models (Dopita & Sutherland 
1996), evidence of ionizing shocks has been found in an in- 
creasing number of AGN from the analysis of their line ra- 
tios (Sutherland et al. 1993; Clark et al. 1998; Koekemoer 
& Bicknell 1998; Villar-Martin et al. 1999c). 

Despite their success, however, shock models remain 
on average less efficient than photoionization models in ex- 
plaining the emission line properties of AGN, especially in 
the UV (Villar-Martin et al. 1997; De Breuck et al. 2000). 
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Therefore, shocks should probably not be regarded as an 
alternative to AGN photoionization, but more likely as an 
additional ionization source that may be dominant only 
under some circumstances. In this paper we investigate 
the possibility that shocks induced by the radio ejecta and 
AGN photoionization coexist in QSOs and radio galaxies. 

To test this hypothesis, we systematically compare ob- 
served emission line ratios, both in the UV and in the opti- 
cal, with models combining the two ionization processes in 
varying proportions. The shock models are computed with 
MAPPINGSIII (Dopita & Sutherland 1996) and the pho- 
toionization models with CLOUDY (Ferland 1996). Our 
approach implicitly assumes that shocks and photoioniza- 
tion don't occur simultaneously in individual clouds. This 
is a reasonable approximation in the case of jet-induced 
shocks located near radio hot spots (Clark et al. 1998; 
Villar-Martin et al. 1999b). At such distances from the 
nucleus, the AGN radiation is weak due to geometrical 
dilution, and the ionization by shocks should be domi- 
nant. Conversely, clouds lying closer to the AGN should 
be mainly photoionized. However, our models probably 
don't apply in the case of shocks induced by inflows or 
outflows (see Contini & Aldrovandi 1983; Contini & Al- 
drovandi 1986; Aldrovandi & Contini 1984, 1985; Contini 
& Viegas-Aldrovandi 1987; Viegas- Aldrovandi & Contini 
1989). 

We also discuss in detail the relationship between the 
radio size and the role of shocks in the ionizing process. 
The present work confirms the conclusion of Best et al. 
(2000b) that large radio sources are dominated by AGN 
photoionization. Here we show that compact sources are 
also dominated by AGN photoionization. Shocks play a 
major role only in some of the intermediate-sized radio 
sources. 

Sect. 2 presents our modeling of the ionizing pro- 
cesses with the help of CLOUDY and MAPPINGSIII. The 
data set used for the comparison with the model predic- 
tions is described in Sect. 3. The results of this compar- 
ison are presented in Sect. 4, with a particular emphasis 
on the NVA1240/HcIIA1640 vs. NVA1240/CIVA1549 dia- 
gram, because of its importance as a metallicity diagnos- 
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tic. In Sect. 5, we discuss the evidence for ionizing shocks 
in radio sources, and investigate the possible origin of the 
link between the shock-photoionization balance and the 
radio size. We summarize our conclusions in Sect. 6. 

2. Modeling the ionizing processes 

2.1. Photoionization by the central AGN 

AGN photoionization models (hereafter A) are computed 
with CLOUDY (Ferland 1996). The emission line spectra 
of photoionized materials depend on the slope of the ion- 
izing continuum, on the ionization parameter £/|], on the 
gas density and on its metallicity. The ionizing continuum 
follows a power law ($„ oc v a ) truncated at 0.01 and 20 

Ryd. 

Two values of the spectral index are considered: a = 
-1, in agreement with the analysis of UV line ratios in 
high- 2 radio galaxies (Villar-Martm et al. 1997), and a = 
-1.5, as suggested by optical line ratios in low-z sources 
(Robinson et al. 1987). The ionization parameter varies 
between log U = -4 and log U = -1. We consider two values 
for the hydrogen density: tih = 10 2 cm -3 and nn = 10 6 
cm -3 . The first value is typical of extended emission line 
regions (McCarthy et al. 1990). The secund one accounts 
for a possible contamination of UV lines by denser regions, 
as suggested by Villar-Martm et al. (1999b). Two different 
metallicities, Z — Zq and Z = 1.5 Zq, are considered. 
The abundances are as in Moy et al. (2001), except for 
nitrogen: in the present work the evolution of the N/O 
ratio with Z follows the empirical relation proposed by 
Bresolin et al. (1999). 

2.2. Ionizing shock 

We use MAPPINGSIII (Dopita & Sutherland 1996) to 
compute the emission line spectra emitted by shocked ma- 
terials. The precursor density is 1 cm~ 3 , as in Dopita & 
Sutherland (1995). We considered shock velocities v be- 
tween 100 and 1000 km s _1 . The magnetic parameter 
B/ \f{n) is fixed to 3 in all calculations. This value cor- 
responds to equipartition between thermal and magnetic 
pressures (Dopita & Sutherland 1996). For the sake of 
consistency, metallicities (Zq andl.5 Zq) and chemical 
abundances are exactly those adopted for photoionization 
calculations. We consider both pure shock (hereafter S) 
and shock+precursor (hereafter SP) models. 

The upward spectrum produced by a v = 1000 km s 
shock is too hard to be fully absorbed by the precursor 
region. To construct an appropriate SP model for v — 1000 
km s , the precursor model is therefore truncated at Phii 
~ 0.7, where F H n is the fraction of ionized hydrogen, so 
that it is "matter-bounded" (Viegas & Prieto 1992). 

1 denned as U = F H ° /n H c, where F H ° is the ionizing flux 
and bh the hydrogen density. 



2.3. Combining photoionization and shock models 

A two-component model is obtained by combining the re- 
sults of a shock model (characterized by the shock velocity 
v) and of a photoionization model (characterized by the 
ionization parameter U) . Only models with similar metal- 
licities are combined in this way. The total flux Fi in a 
given line i is computed as follows: 

Fi = A A/S x F? + (1 - A A/S ) x if A A/s e[0,l] (1) 

where F? and F± are the fluxes predicted by the shock 
and photoionization models, respectively. By varying the 
Aa/s parameter between and 1, we obtain full sequences 
between the "pure shock" and the "pure photoionization" 
cases. 

3. Data sources 

Our data compilation, listed in the Appendix (Table A.l), 
includes emission lines for 369 radio galaxies. 122 have a 
Fanaroff & Riley (1974) classification (10 FRls and 112 
FR2s). The sample also includes 23 compact steep spec- 
trum (CSS) objects and 59 QSOs. For objects with both 
narrow and broad lines, only the narrow component is con- 
sidered. We classify the observations according to their 
spatial coverage: nuclear (including only the central re- 
gions), extended (including only the off- nuclear regions), 
and spatially integrated (covering both the nuclear and 
extended regions). We classified the observations of Simp- 
son et al. (1996) as "spatially integrated" in spite of their 
limited spatial coverage, since they include at least part 
of the extended emission. The radio sizes are estimated 
from the apparent largest angular sizes (LAS), assuming 
Q = 1) A = and H = 50 km s -1 Mpc" 1 . The source 
papers for the LAS are listed in column 3 of Table Al. For 
sources from Morganti et al. (1993), the LAS were directly 
measured from the radio images. 

4. Diagnostic diagrams with combined models 

Figures 1 to 6 present the comparisons between the model 
predictions and our compilation of line ratios. Since large 
ranges of ionization parameter and shock velocity are con- 
sidered, the number of Aa/s sequences is very high. For 
the sake of clarity, only a few sequences, characterized 
by the parameter pair (logC/, v), are plotted on each 
figure. These are (logU» = (-1,300), (-1,100), and (- 
1,1000) for A A / S sequences combining one photoioniza- 
tion model with one pure shock model, hereafter A+S 
sequences (thick solid lines). The first of these sequences 
provides a good fit to the data in many optical diagrams 
and is therefore regarded as the "reference" sequence. The 
latter two sequences usually define the borders of the area 
covered by two-component models. 

One example sequence which includes the emis- 
sion from a precursor (hereafter A+SP sequence), with 
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Fig. 1. Comparisond between the predictions of two component A+S models (thick solid lines) and observational 
data in the three diagrams of Veilleux & Osterbrock (1987). The metallicity is solar, and the density tir- is equal 
to 100 cm~ 3 . The values of log/7 and v are indicated on the plot. Two photoionization model sequences (-4 < 
\ogU < -1), with spectral index a = -1 (thin solid line) and a = -1.5 (dotted line), are shown for comparison. 
An example of an A+SP model sequence with logU ~ -1 and v = 300 km s _1 is also plotted (dashed line). 
Filled squares indicate emission line ratios for spatially integrated (nuclear+extended) emission line regions, 
filled triangles for nuclear regions only, and open circles for extended regions only. Note the similarity between 
the A+S sequences and the pure photoionization sequences. 



(\ogU,v) — (-1,300), is also plotted (dashed line). A "clas- 
sical" AGN photoionization model sequence with spectral 
index a = -1 and -4 < \ogU < -1 is shown for comparison 
(thin line). If not specified, the metallicity is always solar, 
and the hydrogen density nn is always 100 cm~ 3 . Some 
diagrams also plot an additional model sequence (dotted 
line) with different parameters, whose values depend on 
the line ratio under consideration (see below). 

4-1- Optical line ratios 

4.1.1. Veilleux and Osterbrock (1987) diagrams 

Figure 1 shows the comparisons between the var- 
ious models and our compilation of data in the 
three diagnostic diagrams of Veilleux & Osterbrock 
(1987): [OIII]A5007/H/3 vs. [OI]A6300/Ha (Fig. la), 
[OIII]A5007/H/3 vs. [SII]A6720/Ha (Fig. lb) and 



[OIII]A5007/H/3 vs. [Nil] A6584/Ha (Fig. lc). On each 
graph, the A+S sequences (thick solid lines) are very 
similar to the pure AGN photoionization sequence (thin 
line). This result demonstrates that a variation of the 
balance between shocks and AGN photoionization mimics 
a variation of the ionization parameter U , which has long 
been the traditional way of interpreting optical line ratio 
diagrams. 

The a = -1 photoionization models (thin solid line) under- 
predict the observed [OI]A6300/Ha, [SII]A6720/Ha and 
[NII]A6584/Ha ratios, sometimes by as much as 0.5 dex. 
In the case of a = -1.5 (dotted line), the situation is even 
worse, due to the lack of high-energy photons and the sub- 
sequent reduction of the "partially ionized zone" . In order 
to reconcile pure photoionization models with the data, 
one could in principle assume either a lower metallicity 
(see Ferland & Netzer 1983) or a higher density. How- 
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Fig. 2. [OIII]A5007/H/3and [OI]A6300/[OIII]A5007 vs. [OIII]A4363/[OIII]A5007. Symbols as 
in Fig. 1. The photoionization sequence with a = -1.5 is not plotted. Remaining curves as in 
Fig. 1: thick solid lines are A+S models; thin solid line are pure photoionization models (-4 
< \ogU < -1), with spectral index a — -1; dashed line are A+SP models. The data trends 
and dispersions are compatible with A+S models. The A+SP sequence (dashed line) agrees 
almost exactly with the v = 1000 km s _1 A+S sequence. Pure photoionization models are 
clearly excluded. 
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Fig. 3. [OI]A6300 / [OIII]A5007 vs. [OII]A3727/ 
[OIIIJA5007. Symbols as in Fig. 1, and curves as in 
Fig. 2. The A+SP (dashed line) sequence is super- 
imposed to the v = 300 km s _1 A+S sequence. As in 
Fig. 1, the A+S sequences (thick lines) are very sim- 
ilar to the pure photoionization sequence (thin line), 
but with a larger range of [OI]A6300/[OIII]A5007. 



ever, the first hypothesis does not solve the discrepancy for 
[Nil] A6584/Ha, as already emphasized by Ferland & Net- 
zcr (1983) and Veilleux & Osterbrock (1987). The second 
hypothesis is ruled out by the high [SII]A6720/Ha ratios 
observed in our sample. Due to collisional de-excitation, 
the sulphur lines would be weaker than observed if densi- 
ties were higher than tir = 100 cm -3 . A+S models (com- 
bining AGN photoionization and pure shocks; thick lines) 
account for most of the data. A+SP models (combin- 



ing AGN photoionization and shocks+precursors; dashed 
line) are acceptable as well. The sequence with (log U,v) 
= (-1,300) provides the best fit to the data on the three 
diagrams. The dispersion around this sequence, especially 
in [OI]A6300/Ha (Fig. la) and [SII]A6720/Ha (Fig. lb), 
requires that the shock velocity varies between 100 and 
1000 km s _1 . 

Most data points lie above log([OIII]A5007/H/3) 
= 0.5, and only five show a low ionization level 
(log([OIII]A5007/H/3) - 0.). These objects (or regions) are 
probably completely dominated by shocks. The fact that 
there is no nuclear region in this group supports our initial 
hypothesis (Sect. 1) that shocks mainly occur far from the 
nucleus. 

4.1.2. The [OIII]A4363/[OIII]A5007 ratio 

The so-called "temperature problem" has long been 
an outstanding question in the modeling of emis- 
sion lines in AGN. Photoionization models predict too 
low [OIII]A4363/[OIII]A5007 ratios (i.e., too low tem- 
peratures) compared to the observations, (e.g. Tad- 
hunter et al. 1989). Dopita & Sutherland (1995) 
showed that pure shock models are able to account for 
log([OIII]A4363/[OIII]A5007) as high as -1. However, these 
models undcrprcdict the ionization level of the gas and are 
unable to account for the [OIIIJA5007/H/3 ratios in AGN 
(see Fig. 7 in Dopita & Sutherland 1995). 

Figure 2 presents the temperature-sensitive 
[OIII]A4363/[OIII]A5007 ratio vs. the two ionization level- 
sensitive [OIII]A5007/H/3 and [OI]A6300/[OIII]A5007 
ratios (Fig. 2a and 2b, respectively). From these plots 
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Fig. 4. [OIII]A5007/H/3 and [OI]A6300/Ha vs. HcIIA4686/H/3. Symbols as in Fig. 1, and 
curves as in Fig. 2. (except for the dotted line; see below). Many extended regions (open 
circles) show high HcIIA4686/H/3 ratios. They can be accounted for by A+SP models if the 
precursor is matter bounded (dotted line). 



it is clear that A+S sequences solve the temperature 
problem. The bulk of the [OIII]A4363 emission originates 
in shocks, while the [OIIIJA5007 line is emitted mainly 
by the photoionized component. An alternative to the 
A+S solution would be to increase the density, but we 
have shown in the previous section that a significant 
contamination of the optical lines by a high-density zone 
is unlikely. 

The data distribution in the [OI]A6300/[OIII]A5007 
vs. [OIII]A4363/[OIII]A5007 (Fig. 2b) is striking. In the 
classical AGN photoionization hypothesis, one would 
expect an inverse correlation between these two ra- 
tios, since [OIIIJA4363 weakens faster than [OIII]A5007 
at low U (thin line). On the contrary, we find that 
[OI]A6300/[OIII]A5007 and [OIII]A4363/[OIII]A5007 are 
correlated. A Spearman's rank test gives a correlation co- 
efficient of 0.41 with a significance level of ~ 5 10~ 7 . This 
means that objects showing the lowest ionization level also 
have the highest electronic temperature. As illustrated in 
Fig. 2b, this trend is qualitatively reproduced by A+S 
models. Again, this result strongly suggests that the obser- 
vational sequence from high-ionization to low-ionization 
radio sources is actually a sequence from photoionization- 
dominated to shock-dominated objects. 

4.1.3. [OI]A6300/[OIII]A5007 vs. [OII]A3727/[OIII]A5007 

This diagram involves two line ratios which are sensitive 
to the ionization level of the gas. They are frequently used 
in the literature as a diagnostic of the ionization source of 
the gas (e.g. Baldwin, Phillips and Terlevich 1981). Again, 
Fig. 3 illustrates: 

— that a variation of the balance between shocks and 
AGN photoionization quantitatively mimics a varia- 



tion of the ionization parameter U in classical line ratio 
diagrams; and 

— that A+S and A+SP models are better able than pure 
photoionization models to account for the data. 

It is clear from this plot that there is a partial degener- 
acy between the ionization parameter of the photoionized 
component and the value of A^/Sy & n d that log?/ could 
be less than -1 in some objects. Even so, shocks are still 
needed to account for the [OI]A6300//[OIII] A5007. 

Interestingly, shock models predict a relatively small 
range of log([OII]A3727/[OIII]A5007) whichis systemati- 
cally higher than the predictions of AGN photoioniza- 
tion models, in contrast to other line ratios highly influ- 
enced by shocks like [OI]A6300/Ha, [NII]A6584/Ha and 
[SII]A6720/Ha (see Fig. 1). The [OII]A3727/[OIII]A5007 
ratio is therefore an inambiguous indicator of the pres- 
ence of ionizing shocks inside a given object. 

4.1.4. The HcIIA4686/H/3 ratio 

It is well known that the HeIIA4686/H/3 ratio cannot be 
accounted for by classical photoionization models. Figure 
4 shows that A+S and A+SP models (thick lines) ac- 
count for most of the data data in the [OIIIJA5007/H/3 
vs. HeIIA4686/H/3 (Fig. 4a) and [OI]A6300/Ha vs. 
HeIIA4686/H/3 (Fig. 4b). However, these models still un- 
dcrpredict HeIIA4686/H/3 by ~ 0.5 dex for a few data 
(open circles), especially on Fig. 4a. 

Binette et al. (1996) showed that photoionized matter- 
bounded (MB) clouds were one solution which could ac- 
count for the highest HeIIA4686/H/3 ratios observed in 
AGN. In this picture, however, MB clouds are necessar- 
ily located near the nucleus. This scenario is unlikely to 
work for our data since most of the problematic points 
correspond to "extended" regions. As an alternative, we 
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Fig. 5. Comparisons between the predictions of A+S and A+SP models and observational 
data in UV line ratios. Symbols as in Fig. 1. Thick solid curves, thin solid curves, and 
dashed curve are as in Fig. 2, except for Fig. 5c (see below). The dotted curve on Fig. 5a 
and 5b shows the A+S sequence with log U = -3 v = 100, if the density of the photoionized 
component is njj = 10 6 cm~ 3 . The contamination of the emission line spectrum by the NLR 
is strongly suggested by these diagrams. In Fig. 5c, logU — -1 has been replaced by \ogU 
= -2 in the A+S sequence. In Fig. 5d, the dotted line shows the A+S sequence with log U 
= -1 and v = 300 km s^ 1 for Z = 1.5 Z c , r 



argue that MB precursors may well explain the existence 
of off- nuclear regions with high HeIIA4686/H/3 ratios. Pre- 
cursors are likely to lie far from the nucleus if, as our hy- 
pothesis predicts, shocks occur mainly near the radio hot 
spots. Moreover, their ionization level is high, thanks to 
the proximity of their ionization source (Dopita & Suther- 
land 1996). Fig. 4 plots the prediction of the v = 1000 
km s _1 SP model with a truncated precursor (see Sect. 
2.2). We conclude that MB precursors, the existence of 
which is discussed in Sect. 5.1, are able to account for 
log(HeIIA4686/H/3) ~ -0.3 in extended regions. 

4-2. UV line ratios 

We have shown in the previous section that the coexis- 
tence of shocks and photoionization in AGN is strongly fa- 



vored by analysis of optical line ratios. The distribution of 
the data in many diagrams is incompatible with pure pho- 
toionization models, while the A+S sequence with log U= 
<~ -1 and v = 300 km s _1 provides a good fit to the obser- 
vations. The scattering of the data around this sequence is 
well accounted for by A+S sequences with 100 km s _1 < 
v < 1000 km s _1 . This result needs to be confirmed in the 
UV. However, such a study is difficult due to the possible 
contamination of the UV lines by the narrow line region 
(NLR), wich has densities of up to nn = 10 6 cm -3 (Villar- 
Martm et al. 1999b). In what follows we show that the UV 
line fluxes observed in radio sources are compatible with 
the coexistence of AGN photoionization and shocks, if the 
contamination by the NLR is taken into account. 

Figure 5 presents the predictions of the A+S 
and A+SP models compared to the UV data. An 
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A+S sequence with a high density (tih = 10 6 
cm~ 3 ) photoionization component (logt/ = -3) is 
also shown for comparison (dotted line) on the dia- 
grams involving density-sensitive line ratios. The data 
distributions are completely accounted for in the 
CIII]A1909/HeIIA1640 vs. CIVA1549/HeIIA1640 (Fig. 5a) 
and CIII]A1909/HeIIA1640 vs. CIVA1549/CIII]A1909 dia- 
grams (Fig. 5b) only if contamination by a dense region 
is assumed. 

We find, in good agreement with De Breuck et 
al. (2000), that the CIII]A1909/[CII] A2326 ratio is a 
good discriminant between shocks and AGN photoion- 
ization. This ratio is plotted in Fig. 5c as a function of 
CIVA1549/[CII]A2326. Once more, two-component mod- 
els appear to be the only way to account for the ob- 
servations. The fin = 10 6 cm~ 3 model is not plot- 
ted on Fig. 5c, since neither CIVA1549/[CII] A2326 nor 
CIII]A1909/[CII]A2326 is sensitive to density. 

Figure 5d shows the results of two-component mod- 
els for NVA1240/HeIIA1640 vs. NVA1240/CIVA1549. This 
diagram is often used as a metallicity diagnostic in 
AGN environments (Hamann & Ferland 1993). The high 
NVA1240/HcIIA1640 and NVA1240/CIVA1549 ratios in 
radio galaxies and quasars are frequently interpreted as 
evidence of largely super-solar metallicities, not only in 
broad line regions (Hamann & Ferland 1999), but also on 
much larger scales (Villar-Martm et al. 1999b, De Breuck 
et al. 2000). Fig. 5d illustrates that extreme metallicity 
effects are not necessarily needed to explain the intensity 
of the NV line. Many data points can be accounted for by 
A+S and A+SP models at solar metallicity. The objects 
showing the most extreme NVA1240/CIVA1549 ratios can 
easily be accounted for by slightly increasing the metallic- 
ity (Z = 1.5 Z e ) in A+S models (dotted line). 

5. Discussion 

5.1. Evidence for ionizing shocks in extended regions 

An extensive analysis of UV and optical line ratios in ex- 
tragalactic radio sources has been presented in the previ- 
ous sections. The results strongly support the hypothesis 
that the ionized gas in radio galaxies and QSOs is excited 
partly by the radiation from the central source and partly 
by shocks. First, the model sequences obtained by varying 
the shock contribution to the H/3 flux (Aa/s) between 
and 100 % provide better fits to the data than classical U 
sequences in all line ratio diagrams. Moreover, the obser- 
vational sequence from high-ionization to low-ionization 
radio sources is naturally accounted for by varying Aa/s- 
The shock velocity lies between 100 km s _1 and 1000 km 
s _1 . The ionization parameter \ogU of the photoionized 
clouds is typically high, ~ -1. 

Our basic hypothesis (Sect. 1) was that ionizing shocks lie 
far from the nucleus inside radio sources. This assumption 
is a posteriori justified by the high [OIIIJA5007/H/3 ratios 
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Fig. 6. Evolution of the ionization-level sensitive 
[OII]/[OIII] ratio as a function of the radio size in 
radio galaxies. Filled squares represent FRII sources 
and open squares FPJ sources while stars represent 
CSS sources. Triangles are the radio galaxies with- 
out known classification. The horizontal lines show 
the ranges of variation of the [OII]/[OIII] ratio in 
the cases of pure shocks and pure photoionization by 
the AGN. An interpretation of this diagram is given 
in Sect. 5.2 and 5.3. 

observed in all the central regions (Fig. 1), and by the 
fact that evidence for shocks is found mainly in extended 
regions. Another indirect indication of the prevalence of 
shocks far from the nucleus are the large HcIIA4686/H/3 
ratios observed in the external parts of radio galaxies (Fig. 
4). These ratios can be accounted for only by models in- 
cluding a truncated (or "matter-bounded" ) precursor. The 
existence of high-£7 MB clouds photoionized by the AGN 
(Binette et al. 1996) is unlikely so far from the central 
engine. In contrast, the hypothesis of MB precursors is 
reasonable. The ionizing photon source (the cooling post- 
shock medium) is very close. At large distances from the 
nucleus, it is possible that there is not enough material to 
absorb all the ionizing photons emitted by the post-shock 
zone. Finally, the hardness of the upward shock spectrum 
itself does not favor 100 % absorption of the ionizing pho- 
tons. 

5.2. The radio size-ionizing mechanism relation 

From Fig. 1 to 5, it is clear that the balance between 
shocks and AGN photoionization differs from one object 
to another. Recently, Best et al. (2000b) presented con- 
vincing evidence that radio galaxies with radio sizes (D) 
larger than ~ 150 kpc are dominated by AGN photoion- 
ization, and smaller sources by shocks. These results are 
based on the decrease of the [CII]A2326/CIII]A1909 ra- 
tio at large radio size, which has been confirmed by De 
Breuck et al. (2000). A confirmation of the Best et al. 
(2000b) conclusions is clearly needed in the optical. The 
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ionization- level sensitive [Oil] A3727/ [OIII] A5007 ratio is 
ideal for such an analysis. The [OII]A3727/[OIII]A5007 
predictions are higher for shock models than for photoion- 
ization models (Fig. 3). In addition, [Oil] and [OIII] fluxes 
can easily be found from the literature, for objects cover- 
ing a wide range of radio sizes, permitting the construc- 
tion of an ionization level-radio size diagram covering four 
decades in D. We plot on Fig. 6 the comparison between 
[Oil] /[OIII] and radio size for our sample. For comparison 
with Best et al. (2000b) and De Breuck et al. (2000), as 
well as for the sake of homogeneity, only the integrated 
data for radio galaxies, both FRI/II and CSS, are plot- 
ted. Note that wc have [Oil] and [OIII] fluxes for 9 CSSs 
(compared to the 23 CSSs included in our sample). The 
so-called "weak line radio galaxies" (e.g. Tadhuntcr et al. 
1998) are not plotted, due to uncertainty about the im- 
pact of the AGN inside these objects. Fig. 6 also shows the 
ranges of [Oil] /[OIII] covered by photoionization models 
with U > 10~ 3 , and by shock or shock+precursor models 
with 100 km s" 1 < v < 1000 km s" 1 . 

Figure 6 confirms that there is a strong relation be- 
tween the radio size and the ionization level of the gas. The 
most striking feature on Fig. 6 is the trend defined by CSS 
sources, namely a correlation between [OII]/[OIII] and D, 
although an inverse correlation is expected if shocks dom- 
inate in the most compact sources. Large reddening ef- 
fects may explain this relation if the degree of reddening 
is a function of radio size. Such a link is very unlikely, 
as demonstrated by Best et al. (2000b). If CSS sources 
are frustrated by a very dense medium (van Breugel et 
al. 1984), collisional de-excitation effects could explain 
the low [Oil] /[OIII] in the CSS galaxies (since the crit- 
ical density is lower for [Oil] than for [OIII]. In princi- 
ple, this could also explain the correlation between D and 
[Oil] /[OIII] if the most compact sources are confined by 
the densest media. 

To date, however, there is more support for the 
hypothesis that CSS are young rather than frustrated 
(dc Vries et al. 1998; Owsianik & Conway 1998; Fanti 
et al. 2000; Snellen et al. 2000). If it is the case, Fig. 6 
strongly suggests that small radio sources are dominated 
by photoionization, then by shocks when they have grown 
up to ~ 1 kpc, and finally by photoionization when D 
> ~ 150 kpc. One exception to this scheme is 3C236. 
This galaxy, the largest FRII known (D ~ 6 Mpc), has 
a high [Oil] /[OIII] ratio. This is possibly due to the 
existence inside 3C236 of an inner, compact (2 kpc) radio 
source (O'Dea et al. 2001). This source, presumably 
shock-dominated, may dominate the integrated emission 
of 3C236. 
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5.3. Why are large and compact sources 
photoionization- dominated ? 

We now investigate the origin of the relation between ra- 
dio size and shock-photoionization balance. From Fig. 6, 
compact /young sources seem to be dominated by pho- 
toionization. The shock contribution then increases with 
radio size. This may be related to the increase of the co- 
coon surface (and hence of the shock area) as the hot spot 
advances, while the photoionization component remains 
constant. Best et al. (2000b) interpret the prevalence of 
photoionization at larger radio sizes as a consequence of 
the shock front having passed beyond the material sur- 
rounding the AGN. Since most sources with D > 200 kpc 
are photoionization-dominated, this would imply that the 
warm gas around AGN typically expands to a radius of <~ 
100 kpc from the nucleus. 

These hypotheses about the increasing importance of 
ionizing shocks with time in small-scale radio sources, and 
their decreasing importance at larger scales, are in very 
good agreement with the evolution scenario of radio emis- 
sion in Gigaherz Peaked Spectrum and FRIIs recently pro- 
posed by Snellen et al. (2000). 

Let us suppose that the ISM density follows n — 
no(xh/xo)~ S , with xo = 1 kpc and no = 100 cm~ 3 . The 
values of <5 are typically <~ 1.5 (Blundcll et al. 1999). The 
H/3 luminosity emitted by the shocks at the radio cocoon 
interface (Bickncll et al. 1997) depends on v 3 = w/1000 
km s _1 , and on the shock 

L% = 1.91 x 10~ 3 vl A1 n A sh (2) 

The shock velocity v depends on the hot spot-nucleus sep- 
aration Xh as: 

/x h \ {S - 2)/:i 

• — (I) (3 » 

where vq is a function of the jet energy flux Fe and of the 
density no following vq oc (F^/no) 1 ^ 3 . After elimination of 
A s h, Bicknell et al. (1997) show that the shock luminosity 

L H0 is s iven b y : 

/ a \ 0.8 / x — (<5— 2)/5 

10" (sh) ^ 8 »°- 2 (S) 

(4) 

With S ~ 1.5, the luminosity Lf] 1 ^ increases only slowly 
as the hot spot advances. In this case the expansion of 
the radio cocoon alone is insufficient to account for the 
[OII]/[OIII]-radio size relation observed in CSSs. To solve 
this problem, the photoionization contribution may pro- 
gressively decrease during the expansion, for example due 
to the sweeping-out of material by the radio jet. 
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6. Summary and conclusions 

In this paper we have presented strong evidence that both 
shocks and AGN photoionization contribute to the ion- 
ization of the extended gas in radio sources. In "classical" 
optical line ratio diagrams, a variation between and 100 
% of the shock contribution to the line emission mimics a 
variation of the ionization parameter of the gas, but the 
fit is better than with pure photoionization models. For 
[OI]A6300/[OIII]A5007 vs. [OIII]A4363/[OIII]A5007, 
[OIIL1A5007/H/3 vs. [OIII]A4363/[OIII]A5007 and 
CIII]A1909/[CII]A2326 vs. CIVA1549/[CII]A2326, the 
data can be reproduced only with a mixture of shocks 
and AGN photoionization. The [OII]/[OIII] vs. radio 
size diagram shows that the most compact sources (D 
< 2 kpc) and the sources larger than D ~ 150 kpc are 
dominated by AGN photoionization. Shocks dominate 
the ionizing process at intermediate radio sizes. This 
extends to small sizes the results of Best et al. (2000b) 
and De Breuck et al. (2001) obtained from UV line ratios. 

These results are in good agreement with the hypothe- 
sis that the radio cocoon and hot spots lead ionizing shocks 
into the ambient ISM. At the beginning of the expansion 
phase, the shock working surface is small and AGN pho- 
toionization prevails. Then shocks become progressively 
important as the working surface increases and/or the 
amount of material swept out by the jet reduces the im- 
portance of the photoionizcd component. Eventually the 
shock front passed beyond the edge of the ISM surround- 
ing the nucleus and AGN photoionization prevails in later 
phases. 

High spatial resolution data are required to identify 
the shock-dominated zones inside radio sources and con- 
firm our results. Integral field spectrographs are clearly the 
ideal instruments for this purpose, mainly when mounted 
on 10-m class telescopes like the VLT. 
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Reference 1 


Objects types 2 


Ref. for radio sizes 3 


Spatial covering 


Fosbury et al. (1987) 


CSS 


Tadhunter et al. (1994b) 


Integrated 


Robinson et al. (1987) 


RG 


none 


Nucleus & Extended 


Saunders et al. (1989) 


RG 


3CR atlas b 


Integrated 


McCarthy et al. (1990) 


RG 


a 


Integrated 


Gelderman & Whittle (1994) 


RG, CSS, QSO 


a 


Integrated 


Tadhunter et al. (1994a) 


RG 


Krichbaum et al. (1998) 


Extended 


Soifer et al. (1995) c 


QSO 


none 


Integrated 


Simpson et al. (1996) 


RG 


none 


Integrated 


Dey et al. (1997) 


RG 


none 


Integrated 


Morganti et al. (1997) 


CSS 


a 


Integrated 


Rottgering et al. (1997) 


RG 


Rottgering et al. (1994) 


Integrated 


Clark et al. (1998) d 


RG 


Baum et al. (1988) 


Integrated 


Koekemoer & Bicknell (1998) e 


RG 


Subrahmanyan et al. (1996) 


Integrated 


Lacy et al. ( 1998a) s 


RG 


none 


Integrated 


Lacy et al. (1998b) 


QSO 


Lawrence et al. (1993) 


Integrated 


Tadhunter et al. (1998) 


RG, CSS, QSO 


Morganti et al. (1993) 9 


Integrated 


Villar-Martm et al. (1998) 


RG 


a 


Integrated 


Best et al. (1999) 


RG 


a 


Integrated 


Lacy et al. (1999) 


RG, QSO 


a 


Integrated 


Murayama et al. (1999) 


QSO 


none 


Integrated 


Simpson et al. (1999) 


RG 


a 


Integrated 


Snellen et al. (1999) 


RG, QSO 


none 


Integrated 


Villar-Martm et al (1999b) 


RG 


Pentericci et al (1999), 
Rottgering et al. (1994) 


Integrated 


Villar-Martm et al. (1999c) 


RG 


Morganti et al. (1993) 


Extended 


Best et al. (2000a) 


RG 


a 


Integrated 


De Breuck et al. (2000) 


RG 


a 


Integrated 


Palma et al. (2000) 


RG 


a 


Integrated 


Schoenmakers (2000) 


RG 


a 


Integrated 


De Breuck et al. (2001) 


RG 


a 


Integrated 



Table A.l. Col. 1: Reference for line fluxes. Col. 2: type(s) of object included in this source (RG: 
radio galaxies; CSS: compact steep spectrum sources; Sey: Seyfert galaxies). Col. 3: Reference(s) for 
radio size. Col. 4: spatial covering, (a): Reference for emission lines and radio sizes are identical, 
(b): 3CR atlas on line: http://www.jb.man.ac.uk/atlas/basic.htm]. (c): Source of the [OIJA6300 
flux in F10214+4724 (Lacy et al., 1998b). (d): Sum of the Nucleus and EELR components, (e): 
Sum of the Center and Tail components, (f): Sum of the "a" and "c" components, (g) Radio sizes 
estimated from the radio maps. 



